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Abstract. In most of the literature on granular gases it is assumed that the resti- 
tution coefficient e, which quantifies the loss of kinetic energy upon a collision is in- 
dependent on the impact velocity. Experiments as well as theoretical investigations 
show, however, that for real materials the restitution coefficient depends signifi- 
cantly on the impact velocity. We consider the diffusion process in a homogeneous 
granular gas, i.e. in a system of dissipatively colliding particles. We show that the 
mean square displacement of the particles changes drastically if we take the impact 
velocity dependence of e into account. Under the oversimplifying assumption of a 
constant coefficient one finds that the particles spread in space logarithmically slow 
with time, whereas realistic particles spread due to a power law. 

1 Introduction 

Granular gases, i.e. gases composed of particles of mesoscopic or macroscopic 
size which suffer inelastic collisions, may be described under certain assump- 
tions by the same quantities as common molecular gases. One can define for 
these systems temperature, density, velocity field, etc., and describe them 
within the framework of hydrodynamics. Due to the inelasticity of the parti- 
cles, however, there emerge new additional time and length scales. 

From a phenomenological point of view, granular gases behave very dif- 
ferently from molecular gases. One observes cluster formation, e.g., [Q-Q and 
vertexes, e.g., during the evolution of an initially uniform granular gas. 
Although plausible explanations of these effects exist, |l]-|^,[| they still lack 
a comprehensive theory. 

Starting from a homogeneous state, all structures in granular gases de- 
velop after some time-lag. Before noticeable inhomogeneities appear the gas 
evolves in the homogeneously cooling state (HCS), when it gradually loses its 
energy according to inelastic collisions. During this regime of the evolution 
the gas is completely described by the current temperature (which decreases 
with time) and by the velocity distribution function. The HCS which precedes 
all further states of evolution, as e.g. clustering, is the most simple state of 
a granular gas. Understanding the physics of the HCS might, therefore, shed 
some light on more complex phenomena such as cluster and vortex formation. 
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In the present study we focus on the diffusion processes in granular gases 
being in the HCS. We consider the case of a constant restitution coefficient, 
which frequently is assumed in literature but turns out to be an oversimplified 
model of a real collision, and impact- velocity dependent restitution coefficient 
as it is valid for viscoelastic particle interaction. 

2 Collision of Particles 

The microscopic dynamics of granular particles is governed by the (normal) 
restitution coefficient e which relates the normal components of the particle 
velocities before and after a collision, = Vi — Vj and v'^ = v' i — v'j by 

. The unit vector e = ry/\r^\ gives the direction of the 



inter-center vector r t j — r< — Tj at the instant of the collision. From the 
conservation of momentum one finds the change of velocity for the colliding 
particles: 

vl = v i --(\+e)(v i ye}e, = Vj + - (1 + e) (v tj ■ e) e . (1) 

For elastic collisions one has e — 1 and for inelastic collisions e decreases with 
increasing degree of inelasticity. 

In literature it is frequently assumed that the restitution coefficient is a 
material constant, e = const. Experiments, e.g. ](|, as well as theoretical in- 
vestigations |Q show, however, that this assumption is not consistent with 
the nature of the inelastic collisions, it does not agree even with a dimension 
analysis ||. The impact velocity dependence of the restitution coefficient 

e ^ eVy ^ has been obtained by generalizing Hertz's contact problem to vis- 
coelastic spheres 0. From the generalized Hertz equation one obtains the 
velocity-dependent restitution coefficient M 



2/5 



with 



2 Y VR' 



T ■ ■ ■ (2) 



(3) 



3m e .ff(l-^ 2 ) ' 

where Y is the Young modulus, v is the Poisson ratio, and A depends on dis- 
sipative parameters of the particle material (for details see 0). The effective 
mass and radius are defined as 

R eff =R 1 R 2 /{Ri+R 2 ) m e ff =m 1 m 2 /(m 1 +m 2 ) (4) 

with R\/2 and mi/2 being the radii and masses of the colliding particles. The 
constants are given by ||, |) 

Cl = -f§l^l i - „ 1.1534 C 2 = -Cf w 0.7982 . (5) 

2 1 /552/5 jr (2i/io) 5 1 w 
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Equation (^|) refers to the case of pure viscoelastic interaction, i.e. when the 
relative velocity |i>r,e| is not too large (to avoid plastic deformation of the 
particles) and is not too small (to allow to neglect surface effects such as 
roughness, adhesion and van der Waals interactions). 



3 Granular Gas Dynamics 

The diffusion equation 

J A (r) = -D A Vn A (r) (6) 

describes the current J A (r) of particles of sort A at position r against the 
density gradient Vn^(r) with D A being the diffusion coefficient. Using the 
continuity equation 

we obtain the "canonical" diffusion equation: 

d -^P=D A VW{r). (8) 

Usually the term diffusion refers to the motion of a species of particles A in 
a "solution" of other particles B. It is also possible, however, that A and B 
are of the same type, distinguishable only by a feature which does not affect 
the mechanical properties, i.e. mechanically A and B are indistinguishable; 
one may think of A and B having different colors. The process of diffusion of 
tagged particles among mechanically identical ones is called self-diffusion. 

The diffusion (and self-diffusion) coefficient is closely related with the 
mean square displacement of tagged particles with time. Assume the tagged 
particles are located at time t = in the origin r = 0. Then one can write 
for the square average displacement of particles at time i: 

([r(t)} 2 ) = j ' dr[r(t)] 2 n A (r,t) j J drn A (r) (9) 

where n A (r,t)/ J drn A (r,t) is the fraction of particles A located at r at 
time t. (Note that J drn A (r,t) = N = const is the total number of tagged 
particles.) Now we multiply both sides of Eq. (||) with r 2 and integrate over 
dr. Using two times integration by parts and definition (||), one obtains for 
3D-systems 

!([r(£)] 2 )=6A4; ([r(tf) = 6D A t . (10) 

Using the kinematic relation r(t) = J v(t')dt', one obtains 

([r(*)] 3 ) = (jT v(t')dt' £ v(t")dt"^ (11) 
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and encounters with the velocity autocorrelation function (VAF) (v(f')v(t")). 
For gases in equilibrium the VAF depends only on the time difference, \t' — 1"\ 
and decays with a characteristic time t v . Using these properties of the VAF, 
one can perform time-integration in ( pd| ) for t 3> t v , i.e. for t — * oo. Taking 
into account ( [lO| ) one obtains the basic relation 

d a = \ j (v(0)v(t))dt. (12) 

Although not stated explicitly, all the above discussion refers to the case 
of equilibrium gases. Granular gases are a priori non-equilibrium systems, 
nevertheless, the concept of the diffusion coefficient may be generalized for 
such systems. Obviously, this refers only to "liquid" or gaseous phases of 
the material |[(| where the particles are highly mobile. As discussed in the 
introduction throughout this article we assume that the granular material is 
homogeneous and isotropic (HCS). 

Another quantity which we will need below is the "granular temperature" 
Tit), which decreases with time due to the loss of kinetic energy according to 
inelastic collisions. If the particles do not lose too much energy in a collision, 
the temperature changes on a time-scale tq, which is much larger than the 
mean collision time r c , i.e. To 3> t c This condition allows for the defini- 
tion of temperature, but imposes some restrictions for the intervals of time 
and values of material parameters. For a full discussion see p| . 

Whereas the diffusion coefficient D for equilibrium systems is just a con- 
stant, the time dependence of temperature causes the diffusion coefficient to 
be time dependent too. Therefore, the natural generalization of the diffusion 
coefficient for nonequilibrium systems is the diffusivity 

^(Z\r(i)) 2 ) = 6 J D{t')dt' . (13) 

The brackets (• • •) denote averaging over the non-equilibrium ensemble, which 
evolution is described by a time-dependent A-particle distribution function 
p(t) (for simplicity we left the same notation as for the equilibrium average). 



4 The Pseudo— Liouville Operator 

To describe the dynamics of the granular gas we use the formalism of the 
pseudo-Liouville operator C |lj,|l3|. This formalism allows to treat the dy- 
namics of a system of particles with hard-core (i.e. singular) interactions 
formally in the same way as if the particles would interact via a smooth po- 
tential. In classical mechanics the time derivative of any dynamical variable 
B reads 

±B({r i ,v i },t) = {B,H} = i£B, (14) 
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where {. . . , 7i} denotes the Poisson brackets, which imply differentiation of 
the Hamiltonian with respect to the coordinates. For singular inter-particle 
potentials, like a hard-core potential, Poisson brackets are not well defined. 
Nevertheless a Liouville-like operator may be defined: 

^ = E^-^-+E% ( 15 ) 

3 1 i<3 

The first sum in ( |l5| ) refers to the free streaming of the particles (the ideal 
part), while the second sum refers to the interactions of pairs of colliding 
particles {i, j} which are described by the binary collision operator fl4|| : 



Tij = a A I d'e [-Vy ■ e j \ Vij ■ e \S [r^ - ae j (6?. - 1] . (16) 

Here a = 2R is the diameter of particles and Vy ■ e = (vi — Vj) ■ e is the 
normal component of the relative velocity of the colliding pair (which, multi- 
plied by the infinitesimal time dt, gives the length of the collision cylinder). 
The Heaviside function 0(x) selects approaching particles and the <5-function 

represents the hard-core interaction. The operator W ■ is defined as 

v 

%f ( r M rj,vu Vj ---)=f (ri,rj,v-, v- • ■ ■) , (17) 

where / is some function of dynamical variables. The after-collision velocities 
of the colliding particles, v[ and ■uj are related to their pre-collision values 
Vi, Vj via Eq. (|l|). 

The pseudo-Liouville operator allows to perform calculations in a very 
elegant way: Formal integration of Eq. ( pT] ) yields for (t > t'): 

B ({n,Vi},t) = e^-^B {{r u vi\,t') . (18) 
With (|l^) the time-correlation function reads 

(B(t')B(t)) = f drp(i / )S(t')e i£(t_i ' ) J B(i / ) , (19) 



where / dr denotes integration over all degrees of freedom and p(t') depends 
on temperature T, particle number density n , etc., which change on a time- 
scale < > t c . In accordance with the molecular chaos assumption at t ~ r c 
the sequence of successive collisions occurs without correlations. If B does 
not depend on the positions of the particles, its time-correlation function 
reads |p| 

(B(t')B(t)) = (B 2 ) t , e-M'IAaCO (t > t') , (20) 

where (• • denotes averaging with the distribution function taken at time 
t'. The relaxation time tb is inverse to the initial slope of the VAF JTa] . It may 
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be found from the time derivative of (B(t')B(t)) taken at t = t'. Eqs. ( fL9| ) 
and ( p0| ) then yield 

-r B \t') = J dr P (t')BiCB/(B 2 ) t , = (21) 

The relaxation time Tg 1 ^') which depends on time via the distribution func- 
tion p(t'), changes on the time-scale t ^> r c . 



5 Velocity Autocorrelation Function 
and Self— Diffusion Coefficient 

5.1 Constant Restitution Coefficient e 

Following the idea described in section || we calculate the velocity autocor- 
relation function and the self-diffusion coefficient. In the previous section we 
discussed a method to calculate the evolution of a general dynamic function 
B(t). Now we specify B(t) to be the velocity of a tagged particle, say vi(t). 
Then with 3T(t) = (v 2 ) t Eqs. @0) (with (|J]l|)) read 

(-ui(i') • i7 X (t)) = 3T(t / )e- |t_t ' l/T " ( *' ) (22) 



r- 1 (0 = (JV-l) \ (23) 

To obtain ( p3|) we take into account that CqVi = 0, TtfVi — (for i ^ 1) and 
the identity of the particles. The calculation of T~ Y it!) may be performed if 
we assume that the distribution function pit') is a product of the coordinate 
part, which corresponds to a uniform and isotropic system, and the velocity 
part being a product of Maxwellian distribution functions 

ni} [27rT(i')] 3 / 2 ' i - i ''--' iV - W 

We want to mention that a more sophisticated analysis (e.g. [p!6j|l7| ) shows 
that the velocity distribution function deviates from the Maxwellian. In this 
paper we neglect these deviations which are small for small inelasticity of the 
particles (e — > 1). Calculation of the diffusion coefficient with non-Maxwellian 
velocity distribution function has been performed in fl8|| . The result (depend- 
ing on the inelasticity) may differ quantitatively from that given here, but 
its functional form remains unaffected. 

Integration over the coordinate part in ( p3|) yields 

[N - 1) j p(t')S fry - <re) dn • • ■ dr N = ng 2 (a) JJ , (25) 
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where we use the definition of the two-particle distribution function and where 

ff2 (a) = i(2-r,)/(l- ?7 ) 3 (26) 

is its contact value (l^,|l5|], which depends on the packing fraction r\ = ^irna 3 . 
With 

[»iTi2«i) = ~ («iaTi2t>i) (27) 



due to the collision rules and definition (16) one finally arrives at 

ng 2 (a)a 2 f f 2 0(-v 12 ■ e)\v 12 ■ e\ (v 12 ■ e) 2 (l + e) 



(28) 

where 

4>{v 12 ) = (4ttT)- 3 / 2 exp(-« 2 2 /4T) (29) 

is the Maxwellian distribution for the relative velocity of two particles. For e 
not depending on v\ 2 Eq. ( p8| ) yields 

r v -\t) = l±l^g 2 (a)^m= e -±^r E 1 (t), (30) 

where Te(t) = | r c (f) is the Enskog relaxation time Jl3|]. For the granular gas 
it depends on time according to the same time-scale as temperature. Accord- 
ing to Eq. ( |30| ) the velocity correlation time for inelastic collisions is larger 
than for elastic collisions. This follows from partial suppression of backscat- 
tering of particles due to inelastic losses in their normal relative motion. As 
a result for inelastic particle interaction the angle between crossing trajec- 
tories after a collision is smaller than before the collision, while for elastic 
interaction both angles are identical. Therefore, the velocity correlation time 
is larger for inelastically colliding particles. 

Using the velocity correlation function one writes 

({Ar(t)) 2 ^=2^ dt'3T(t') £ *"e-l'"-*'l/ T " (t ') . (31) 

On the short-time scale t ~ r c , T(t') and T v (t') may be considered as con- 
stants. Integrating in ( |3l"| ) over t" and equating with ( [l3| ) yields for t ^> r c ~ 
t v the diffusivity (time-dependent self-diffusion coefficient) 

D(t)=T(t) Tv (t). (32) 

Using the pseudo-Liouville operator one can also describe the time-depen- 
dence of the granular temperature: From ( p"4| ) it follows (see also JL9[): 
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Calculations similar to that for T v {t) yield 

"4 



T(t) 



yi 2 g 2 {a)a 2 n^-t 2 ) 



rp3/2 



So that finally one obtains a well known result for temperature 
T(t)/T = [l+ 7o i/T c (0)p 2 , 



(34) 



(35) 



where 70 = (1 — e )/6 |LLp], To is the initial temperature at t = 0, and 



7- c (0) =47r L ^g 2 (a)a z n^/T 



(36) 

is the initial mean collision time. Thus, using Eqs. (^0|) and ( |32"| ) one obtains 
for the self-diffusion coefficient for the case e = const. 



D(t) 



6+1 



D [l+ lo t/T c (0)} 1 



where D is the Enskog self-diffusion coefficient for elastic particles: 

n-l ^1/2 2 ( \ m-1/2 

D = -7T ' a g 2 {<j)nT a 
Correspondingly, the mean-square displacement reads for t 3> t c (0): 

{Ar{t)f)~\ogt. 
5.2 Impact Velocity Dependent Restitution Coefficient 



(37) 

(38) 
(39) 



As mentioned before the assumption of a constant restitution coefficient con- 
tradicts basic physical understanding and does not even agree with a dimen- 
sion analysis |^| . Nevertheless it was used in many publications for simplicity 
of the calculation. In this section we will calculate the diffusivity for a gas of 
viscoelastic particles for which e is given by Eq. (Q) and, surprisingly, we will 
find a qualitatively different result. 

The velocity correlation time t v can be found from Eq. ( |28| ) where e 
depends on the impact velocity as described by Eq. (|^) : 



r-\t) 



3 p /21 
4 V 10 



27 
4lT 



-r 



dAa 2 ' 5 (4T(t)) 1/10 + 
^ C\A 2 a^ {4T{t)) 1/r ° T • • • , 



(40) 



where P (x) is the Gamma-function, te is given by Eq. fl30| ) and we use 
Eq. (||), which relates the coefficients C\ and C 2 . 

Performing calculations similar to those for the constant restitution coef- 
ficient, one obtains for the time evolution of temperature 



T/T - -(5/'S)t^ [(T/T ) 8 / 5 - bS (T/T0) 17 / 10 + 



(41) 
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where we introduce the small (dissipative) parameter 

5 = Aa 2 / 5 T^ /W , (42) 
the characteristic time for the temperature evolution 

r^ 1 = S ■ r^O)- 1 • | 2 1 / 5 C ± r (j^J « 0.831928 ■ S ■ r^O)" 1 (43) 

and the numerical constant b = (18/5)2 1 / 5 Ci J T(l/5)/r(l/10) w 2.3018. Solv- 
ing Eq. ( ftl| ) and expanding the result in terms of the small parameter 8 yields 

T(t) = (l + -) 3 + ai s(l + -) 6 +a 2 S 2 + 2 + ■■■ (44) 



To \ T J \ TQ J \ TO 

with 01 and ai being pure numbers J20I . 

From Eqs. (^ , ^4] , ^0|) follows for the self-diffusion coefficient: 

m = (l + l - Y + a 3 5(l + IV 1 + c**' f 1 + IV* + • • • , (45) 



-Do \ Toy \ Toy \ T o 

where Do is given in Eq. ( |38| ) and 03, 0,4 are pure numbers pcj ]. Correspond- 
ingly, the mean-square displacement reads for tq <C t: 

((Mt)f) ~t 1/6 +a 3 «J log*. (46) 
This dependence holds true for 

.(O)^ 1 «t«r c (0) ( 5- 11 / 5 , (47) 



TV 



where the first inequality follows from the condition tq <C t, while the second 
one follows from the condition T c (t) -C tq which makes the concept of the 
temperature meaningful. 

Comparing Eqs. ( flrf ) and (|3^) one notes that the impact-velocity depen- 
dent restitution coefficient (||) leads to a significant change of the long-time 
behavior of the mean-square displacement of particles in time. Compared 
to its logarithmically weak dependence for e — const, (simplified collision 
model), the impact-velocity dependence of the restitution coefficient (g) as it 
appears for viscoelastic particles gives rise to a considerably faster spreading 
of particles according to a power law. 



6 Results and Discussion 

We studied the diffusion of particles in a homogeneously cooling granular 
gas. With the assumption of molecular chaos we calculated the velocity time- 
correlation function and the time dependent self-diffusion coefficient (diffu- 
sivity) D(t). For the case of constant coefficient of restitution e = const, the 



116 N. V. Brilliantov and T. Poschel 



diffusivity D{t) is expressed in terms of the model parameter e. For the more 
realistic case of an impact velocity dependent restitution coefficient, as it has 
been derived for viscoelastic particle interaction, we found a relation which 
expresses the diffusivity in terms of material constants of the particles and 
characteristics of the granular gas, such as temperature, density, etc. 

For granular particles suffering viscoelastic collisions we found that the 
mean-square displacement grows with time as a power law ((Ar(t)f^ ~ i 1 / 6 , 

i.e. much faster than the logarithmic growth ^(Z\r(i)) 2 ^ ~ logi for the case 

of a constant restitution coefficient. It worth to note that qualitatively this 
power-law dependence (as well as the logarithmic one) simply follows from 
scaling arguments and the time-dependence of temperature. Indeed, the av- 
erage velocity scales as v ~ T 1 / 2 , and therefore as ~ t~ x for the constant 
restitution coefficient and as ~ i -5 / 6 for the impact-velocity dependent coef- 
ficient. The diffusivity in granular gas scales as D ~ 1 2 /t c , where I ~ <r -2 n -1 
is the mean-free path, which does not change with time (in the regime preced- 
ing clustering), and r c ~ l/v is the mean-collision time. Thus, D ~ Iv ~ T 1 / 2 , 
and we obtain that the mean-square displacement, J D(t)dt, scales as ~ logt 
in the former case and as ~ i 1 / 6 in the latter. 

What will be the impact of this apparently dramatic difference in the 
time dependence of (^(Ar(t)) 2 ^ for the properties of granular gases? Obvi- 
ously, viscoelastic particles spread (and, therefore, mix) quicker than particles 
which interact via e = const. Since the temperature decreases more slowly 
for the former case, as ~ i -5 / 3 , as compared with ~ t~ 2 for the latter one, 
retarded clustering may be expected. Let us explain this in more detail: Ac- 
cording to the linear stability analysis pi Q , long-range density fluctuations 
in homogeneous granular gas occur to be unstable, i.e. they grow exponen- 
tially with time, leading to cluster formation. The critical value of the stable 
wavelengths A c = 2n/k c (k is the wave number) reads A c ~ l/y/l — e 2 jjj, 
see also 0. All density fluctuations with A > A c are unstable. For the case 
of viscoelastic particles the effective value of e permanently grows with time 
(as temperature decreases), and A c grows accordingly. This implies that the 
size of the regions, where the gas is still homogeneous, grows with time, since 
density fluctuations do not blow up on length scales smaller than A c . Thus 
for the velocity-dependent restitution coefficient the critical wavelength in- 
creases with time (i.e. one has a time-dependent criterion for the size of the 
region which is unstable in the HCS). As compared to the case e = const., 
clustering occurs on a larger length-scale; hence the conditions to remain in 
the homogeneous medium persist for a longer time for viscoelastic particles. 

Obviously, self-diffusion directly counteracts clustering. Therefore, we con- 
clude, that for a granular gas of viscoelastic particles which implies an impact 
velocity dependent restitution coefficient, clustering is retarded as compared 
with a gas under the oversimplified assumption e = const. 

We thank M. H. Ernst and I. Goldhirsch for valuable discussions. 
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